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Abstract: We show that the naive apphcation of the Kibble mechanism seriously 
underestimates the initial density of cosmic superstrings that can be formed during 
the annihilation of D-branes in the early universe, as in models of brane-antibrane 
inflation. We study the formation of defects in effective field theories of the string 
theory tachyon both analytically, by solving the equation of motion of the tachyon 
field near the core of the defect, and numerically, by evolving the tachyon field on 
a lattice. We find that defects generically form with correlation lengths of order 
rather than H^^. Hence, defects localized in extra dimensions may be formed 
at the end of inflation. This implies that brane-antibrane inflation models where 
inflation is driven by branes which wrap the compact manifold may have problems 
with overclosure by cosmological relics, such as domain walls and monopoles. 
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1. Introduction 

Although it is notoriously difficult to test string theory in the laboratory, excit- 
ing progress has been made on the cosmological front through the realization that 
superstrings could appear as cosmic strings within recent popular scenarios for brane- 
antibrane inflation Gravity wave detectors and pulsar timing measurements 
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could thus give the first positive experimental signals coming from superstrings 
[§] (see, however, p). If seen, it will be challenging to distinguish cosmic super- 
strings from conventional cosmic strings. One distinguishing feature is the possibility 
that superstrings have a smaller intercommutation probability than ordinary cosmic 
strings 0. In this paper we consider whether the mechanism of formation of cos- 
mic superstrings might provide another source of distinction, by studying in detail 
the formation of string and brane defects, emphasizing differences with the standard 
picture of defect formation. 

To set the stage, let us briefly recall how defect formation works in a standard 
theory |Q (see, for example p|, for modern review). Consider a scalar field theory 
with potential f (|0p — ^7^)^- In the standard picture as the universe cools below some 
critical temperature Tc the U{1) symmetry is broken and (f) rolls to the degenerate 
vacua of the theory |0| = 77. Causality implies that the field cannot be correlated 
throughout the space and hence the field rolls to different vacua in different spa- 
tial regions leading to the formation of topological defects. The defect separation 
is set by the correlation length, ^. For a universe expanding with Hubble rate H 
one expects C, < and hence a minimum defect density of about one per Hubble 
volume. However, this is only an upper bound on ^. A more careful estimate can be 
made using condensed matter physics methods: equating the free energy gained by 
symmetry breaking with the gradient energy lost. Very close to Tc thermal fluctua- 
tions which can restore the symmetry are probable; however, once the universe cools 
below the Ginsburg temperature, T^, there is insufficient thermal energy to excite a 
correlation volume into the state = and the defects "freeze out." For the scalar 
field theory described above this estimate yields a correlation length of order the 
microscopic scale: ^ ~ X^^t]"^. 

The physics of tachyon condensation on the unstable Dp brane- ant ibrane system 
is quite different however, due to the peculiar form of the action P| 



The tachyon potential has a "runaway" form and there are no oscillations of the field 
near the true vacuum which can restore the symmetry. The decaying exponential 
potential of the complex tachyon field multiplies the kinetic term. Once condensation 
starts, T rolls quickly to large value, and damps the gradient energy exponentially. 
This essentially eliminates the restoring force which would tend to erase gradients 
within a causal volume in an ordinary field theory. 

In this paper we perform a quantitative analysis of the formation of string de- 
fects starting from the unstable tachyonic condensate that describes unstable brane- 
antibrane systems. Having established that defects form with an energy density 
comparable to that which is available from the condensate, we then examine the 
possible cosmological consequences of this larger-than-expected initial density. 
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The organization of this paper is as follows. In section |^ we briefly discuss brane- 
antibrane inflation and the formation of defects at the end of brane-antibrane infla- 
tion by tachyon condensation. In section |^ we study the formation of codimension- 
one branes in the decay of a nonBPS brane both analytically as well as through 
lattice simulations. In section ^ we repeat the study, this time for the formation of 
CO dimension-two branes in the decay of the brane-antibrane system. We contrast 
these results with the formation of conventional cosmic strings in section ^ The 
reader who is interested in the cosmological implications may wish to skip directly 
to section ^ where we consider new constraints on brane inflation models which may 
arise due to the large initial density of defects. Conclusions are given in section ^ 
In an appendix we justify the assumed initial conditions for the tachyon fluctuations 
which lead to defects. 



2. Brane Inflation 

Relic cosmic superstrings can form at the end of inflation driven by the attractive 



interaction of branes and antibranes |]10[- In this picture, inflation ends when the 



brane and antibrane (or a pair of branes oriented at angles) become sufficiently close 
that one of the open string modes stretching between the branes becomes tachyonic. 
The subsequent rolling of this tachyon field describes the decay of the brane-antibrane 
pair. Quantum fiuctuations produce small inhomogeneities in the tachyon field, 
which will cause it to roll toward different vacua in different spatial regions, leading 
to the formation of topological defects which are known to be consistent descriptions 
of lower dimensional branes. 

The branes which drive inflation must span the three noncompact dimensions 
and may wrap some of the compact dimensions. The defects which form are lower 
dimensional branes whose world-volume is within the world-volume of the parent 
branes. The argument was made in (illustrated in flgure |1|) that by applying 
the reasoning of the Kibble mechanism to the formation of the lower-dimensional 
branes one concludes that the branes which are copiously produced as topological 
defects will always wrap the same compact dimensions that the parent branes wrap 
and hence these branes appear as cosmic strings to the 3-dimensional observer.^ This 
argument is based on the fact that the size of the compact dimensions is orders of 
magnitude smaller than the Hubble distance during (and at the end of) inflation and 
hence there are no causally disconnected regions along the compact dimensions (a 
reasonable estimate is H^^ ~ 6 x lO^AI^^ while typically the extra dimensions are 
compactifled at a scale R closer to M~^). Therefore along these directions the tachyon 
fleld would always roll toward the same vacuum and no topological defects would 



^The production of phenomenologically dangerous monopole-like and domain wall- like defects is 
suppressed. 
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Figure 1: Illustration of the Kibble argument in the case of branes spanning extended 
dimensions and wrapping compact ones. The compact dimension is much smaller in size 
than the Hubble distance, l/i? ^ ^y, so there are no causally disconnected regions along 
this dimension. The branes that will form as topological defects will wrap the same compact 
dimensions as the parent branes, and will be localized on the extended dimensions. 



form. The causally disconnected regions occur only along the extended directions, so 
the defects are localized along the extended directions. An identical argument implies 
that cosmic strings should form with a density of about one string per Hubble volume. 

Reference [U uses the analogue of the condensed matter physics argument para- 
phrased in section ^ to estimate the correlation length of the tachyon field at the 
end of inflation. Since the universe has zero temperature at the end of inflation the 
thermal fluctuations are replaced by the quantum fluctuations of de Sitter space: 
H/{27i). In that analysis the potential depends on the brane separation and as a 
result of the relative motion of the brane and the antibrane, the curvature of the po- 
tential at T = changes sign. It is also assumed that the tachyon potential, V{T), 
has a minimum for some finite value of T. In this case there is an interplay between 
the correlation length, C,c given by the curvature of the potential at T = 0, and the 
Ginsburg length, which is the largest scale on which quantum fluctuations of the 
field can restore the symmetry. As the brane separation decreases and the shape of 
the potential changes, the correlation length and the Ginsburg length change as well. 
The density of defects freezes when = ^g- This analysis yields a correlation length 
which is not much different from the Hubble scale: ^ ~ 1.6 x 10^M~^ ~ 2.7H~^. 

There are several reasons to consider a more quantitative study of defect for- 
mation at the end of brane inflation. Firstly, the estimate of one defect per Hubble 
volume provides only a lower bound on the defect density via causality. The actual 
network of defects produced is determined by the complicated tachyon dynamics 
as it rolls from the unstable maximum of its potential to the degenerate vacua of 
the theory. Second, the effective field theory which describes the dynamics of the 
tachyon ^ has a rather unusual causal structure |11| and the usual reasoning of the 
Kibble mechanism may not be applicable. In [|ll[] it was found that in the case of the 
homogeneous rolling tachyon the small fluctuations of the field propagate according 
to a "tachyon effective metric" which depends on the rolling tachyon background. 
As the tachyon rolls towards the vacuum the tachyon effective metric degenerates, 
the tachyon light cone collapses onto a time-like half line and the tachyon fields at 
different spatial points are causally decoupled. 
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Finally we comment on the validity of the effective field theory description. At 
some point the effective field theory description of the decaying brane will become 
inappropriate and the correct degrees of freedom will be the decay products of the 
brane annihilation. Since the topological defects of the tachyon field form in a short 
time of order M~^, |12[-|I1], the effective field theory description should be applicable 



during the period of defect formation. Furthermore, there are reasons to believe that 
the effective field theory is valid even at late times when the tachyon is close to 
the vacuum, since in this regime the tachyon effective field theory may give a dual 
description of the closed strings which are produced during brane decay O] . 



3. Tachyon Kink Formation in Compact Spaces 



In this section we consider the formation of codimension-one branes from tachyon 
condensation on a nonBPS brane. We study defect formation both by solving the 
full nonlinear equations of motion on a lattice, as well as providing analytical ap- 
proximations to the behavior of the tachyon field in different regimes of interest. We 
study the evolution of the tachyon field starting from a profile which is close to the 
unstable vacuum T = (in this respect our analysis is similar to ||16| , p!7|). In our 
analysis there is no parameter which causes continuous variation of the potential, 
we consider the brane- ant ibrane system to be coincident at t = when the initial 
conditions are imposed (this potentially important difference should be kept in mind 
when comparing our results to [|l|). 



3.1 Action and Equation of Motion 



We will work with the tachyon Dirac-Born-Infeld effective action & 18, 19 



(3.1) 



The action (|3.1| ), with V{T) = Vo / cosh(T/v2a;'), can be derived from string theory 
in some limit [^] and has been shown to reproduce various nontrivial aspects of the 



full string theory dynamics [19, 21]. Here we take the potential to be 



V{T)=T,exp{-Tyb') (3.2) 

where Tp is the Dp-brane tension. The constant b determines the tachyon mass in 
the perturbative vacuum as Mt = a/2&^^ ~ Qualitatively the results will 

depend very little on the specific functional form of V{T). In fact, for much of 
the analytical analysis which follows we will not even need to make reference to the 
specific functional form of V^(T).^ We consider real tachyon fields in a spacetime with 

^Provided of course that V'{T^O)^0 and V{T ^ ±oo) 0. 
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three expanding, noncompact dimensions {x*} (where i = 2,3,4) and one compact 
dimension = x which we take to be static. The metric is 



gf^^dx^'dx" = -df + dx^ + a{tf6ij dx' dx^ . (3.3) 

For simphcity we take the tachyon field to depend only on time, x° = t, and the 
compact spatial coordinate x^ = x. The equation of motion which follows from ([37 
with the metric (|3.3|) is 



(3.4) 



(1 + T'2) f = T" (l - f^^ - (^3Hf + - + T'2 j + 2fT'f' 

where f = d^T = dtT, T' = diT = d^T, V'{T) = dV/dT and H = d/a. 
3.2 Lattice Simulation of Kink Formation 



We have solved (3^) on a lattice for different values of the compactification radius 



and the Hubble parameter H (which we take to be constant for simplicity). We 
consider vanishing initial velocity T{t = 0,x) = Ti{x) = 0. The initial profile 
T{t = 0,x) = Ti{x) is a Fourier series truncated such that the minimum wavelength 
is comparable to the lattice spacing.^ The amplitudes of the Fourier coefficients are 
given by a random Gaussian distribution and the overall amplitude of the initial 
profile is chosen to be small compared to b. 



As in [|T3| we find that the gradient of the tachyon field near the core of the kink 
becomes infinite in a finite time, forcing us to halt our evolution. After this point the 
codimension-one branes have formed and if one wants to follow the evolution beyond 
this time it is necessary to consider the dynamics between branes and antibranes in 
a compact space; the field theory description is no longer adequate. 

Typically the initial profile crosses T = at many points. In the early stages 
of the evolution the field begins to grow due to the small displacement from the 
unstable vacuum. During this phase of the evolution many of the small fiuctuations 
of the field will straighten themselves out. Large Hubble damping tends to kill off 
the high frequency fluctuations faster. At the end of this initial stage there may be 
several locations where the fleld stays pinned at T = 0, depending most crucially on 
the radius of compactiflcation. The evolution very quickly enters a nonlinear regime 
in which the field begins to roll quickly towards T —>■ +oo where T > and T — oo 
where T < 0. 

The most important parameter for determining the density of the defect network 
once the daughter branes have formed is the radius of compactification, R. Perhaps 
surprisingly, we find that for a compactification radius as small as a few times b, 



^See the appendix for a discussion of the vahdity of these initial conditions. 



- 6 - 



tachyon kinks will form. Once the field enters the nonlinear regime the defect for- 
mation depends only on the local physics near the core of the kink. In figures 0, ^ 
we plot the tachyon field versus t and x, showing the formation of co dimension-one 
branes from the decay of a nonBPS brane for radii of compactification R = 8My ^ 
and R = 15M^^. The Hubble constant is taken to be vanishing, H = 0, in these 
figures. We have also considered H ^ and find that the Hubble damping has little 
effect on the final kink/anti-kink network. Hence we find that tachyon kinks do form 
in the compact directions even when the field is initially in causal contact throughout 
the extra dimension. 




Figure 2: Formation of tachyon kinks 
for R = SM^^^ . The time axis is mea- 
sured in units of M^^ . 



Figure 3: Formation of tachyon kinks 
for R = 15M^^ . The time axis is mea- 
sured in units of M^^ . 



3.3 Analytical Study of Tachyon Kink Formation 

Here we describe analytically the formation of tachyon kinks. Since the full equation 
of motion (|3.4|) is difficult to solve analytically for arbitrary initial data, we study 



the dynamics of (|3.4| ) in different regimes: near the core of the defect, where the field 
stays pinned at T = 0, and away from the core, where T rolls towards the vacuum. 
To simplify our analysis we neglect the compact topology of the extra dimension x, 
which should be a good approximation since the defect solutions are infinitely thin 
and therefore highly localized. 

3.3.1 Solutions Near the Core of the Defect 

Here we briefly review the analytical studies of kink formation near the core of 



the defect presented in [0. Consider initial data T{t = 0,x) = Ti{x) and T{t = 
0,x) = Ti{x) = 0. The field should start to roU where Tj(x) ^ due to the small 
displacement from the unstable maximum V'{T) = 0, and furthermore it must stay 
pinned at T = at the core of the kink. At t = the equation of motion ( |3.4| ) is 

T.(x) (1 + T^ixr) = T^ix) + ^r,(x) (1 + T[{xf) . 
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Clearly any point xq where Ti{xo) = T"{xo) = will be a fixed point where T(t, xq) = 
= T{t,xo) throughout the evolution. We restrict ourselves to considering only 
initial data such that sgn(Tj(a;)) = sgn(Tj(x)) for all x to ensure that the solutions 
are increasing. To analytically study the dynamics near xq it is therefore reasonable 
to make the ansatz 

T{t,x)^u{t){x-xo). (3.5) 

For u{t) = const oo this solution corresponds to the singular static soliton solution 
of Sen |TP|. Inserting the ansatz (|3.5|) into (|3.4|) and working only to linear order in 



{x — Xq) one obtains an ordinary differential equation for the slope 

it = + 2 /"^ — 3Hu. (3.6) 
0^ 1 + 

We have solved ( ^.6| ) numerically for both constant H and ~ l/(t + to)- We find 
that generically the solutions of (|3.6|) become singular in a finite time t = tc- Larger 
H tends to delay the onset of the singularity and soften the singular behaviour. To 
understand this finite time slope divergence analytically it is simplest to consider 
H = 0. In this case if we assume that initially ii is close to zero then the second term 
on the right hand side of ( p.6|) can be neglected and therefore 

„(t)^M+e^*/f + „_e--^*/'' 

at early times. Clearly ii grows quickly and the second term on the right hand side 
of ( p. 61) very quickly becomes important. In the regime where u and its derivatives 
become large the solution has the behaviour 

u{t) ^ 



tc-t 

where the critical time tc depends on the initial conditions. We find that within 
finite time the slope becomes singular and the time dependent tachyon field near the 
core of the kink coincides with the singular soliton solution of Sen Hence we 
conclude that the codimension-one brane forms in a finite time. 

This finite-time slope divergence was observed both numerically and analytically 
in and leads to a finite-time divergence in the stress-energy tensor. This diver- 



gence was also found in an exact string theory calculation in [O. For the potential 



2[) the stress-energy tensor near x = Xq and t = t^, with if = 0, is 1 14 



^00 = exp 



tc-t "V b^itc-ty 



Hence, with the normalization b = lyTxa' proposed in |T9|, as t ^ tc we have 
^00 Tp-i5(a; — Xo). Similarly one can show that in this regime Tqi = T\\ = and 
that Tij = —Tp_i6ij6{x — Xq). In the limit of condensation, then, the stress-energy 
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Figure 4: Behavior of the slope at the core of the kink for different initial conditions. The 
time axis is measured in units of b. 



near x = Xq is identical to that of a D(p — l)-brane. If we take into account the 
roUing of the tachyon for x ^ xq then there will be an additional component to T^jy 



corresponding to tachyon matter, as in 22 



3.3.2 Solutions Near the True Vacuum 

Away from the site of the kink the field will roll towards the true vacuum T ±oo 
so that V(T) — > at late times for x xq. As V(T) — > one can show both 
analytically and numerically that ( |3.4] ) has an attractor which satisfies 



d^Tdf.T + 1 



0. 



(3.7) 



This attractor holds for arbitrary metric and hence the solutions we describe below 
hold equally well for vanishing or nonvanishing H. Exact solutions of (p.7|) were found 
in for arbitrary initial data. The solution is given along a set of characteristic 
curves which, in 1 + 1 dimensions, are defined by 



x{q,t) = q 



■.t 



(3.1 



where the parameter q defines the initial position of the curve on the x axis such 
that x{q,t = 0) = q and Tj(x) is the Cauchy data at time t = 0. The parameter q 
should be thought of as labeling the curves. The value of the tachyon field along a 
given characteristic curve is 



T(g,t)=T,(g) + 



t 



(3.9) 



V^ + V{qy' 

For generic initial data these solutions form caustics in finite time where the second 
derivatives of the field blow up and the effective description (|3.1|) breaks down. In 



the present context the caustics lead to singularities in the stress-energy tensor [25 



for X Xq. Caustic formation does not arise in our lattice simulations since the 
finite-time first derivative divergence associated with defect formation forces us to 
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halt our evolution, and this generically seems to occur before caustic formation. The 
formation of caustics is predicted by a wide range of tachyon effective field theories 
[ p^ and is presumably not present in the full string theory dynamics. 



4. Tachyon Vortex Formation in Compact Spaces 



In this section we generalize the results of section § to consider the formation of 
codimension-two branes from tachyon condensation on the brane-antibrane pair. 
This situation is of more direct interest to brane-antibrane inflation since inflation 
ends when a brane-antibrane pair annihilate. This situation is also more realistic 
since the stable branes in a given string theory are those whose dimension differs by 
a multiple of two. 

4.1 Action and Equations of Motion 

We wish to generalize the results of the previous section to consider complex tachyon 
fields which depend on time and two spatial coordinates which we compactify on a 
square torus. We generalize the action (p.l|) to 



S 



-2Tp I e 



1 + d^Tdt'T 



-9 d 



'3+2+1 



X 



(4.1) 



with metric 



g^^dx^dx" = —dt^ + dx^ + dy'^ + a{t)'^5ijdx'^dx^ 

where x^ = x and x"^ = y are Cartesian coordinates on the torus and denotes 
the three noncompact dimensions, as in section ^. We consider tachyon fields which 
depend on x'^ = t, x^ = x and x^ = y. If we separate the tachyon field into into real 
and imaginary components as 

T{t,x,y) = Ti{t,x,y) + iT2{t,x,y) 

then ( [4.1|) may be rewritten in terms of components as 

S = -2Tp J exp {-TjTj/ b^) v^l + d^Tjd^^Tj ^ d^+^+'x 

where the upper-case Roman indices label the real and imaginary components of the 
tachyon field (J = 1, 2). The equation of motion is (no sum over J, sum over J and 
K): 



dVi\T\) 
dTi 



-gg^^dpTi 



+ 



dV{\T\) 



^^/l + d^Tjd-Tj = 0. 



VI + d^Tjd-Tj 



(4.2) 
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We note that the action (|4.1|) has not been derived from first principles and 
has several drawbacks from a theoretical perspective. For other proposals of the 
effective field theory on the brane-antibrane pair see One theoretical 

difficulty is that the action ( [4 .11) does not evade Derrick's theorem |3^. This is 
of no practical concern to us since we study time dependent solutions and since our 
interest is in defect formation in a compact space where charge conservation precludes 
the possibility of isolated defect solutions. More seriously, for the action ( |4.1| ) a 
static profile of the form T = u{x + iy) with u ^ oo does not correctly reproduce 
the stress tensor for a codimension-two brane. We will attempt to address this 
difficulty by considering an alternative description of the complex tachyon dynamics 
in a subsequent subsection. 



4.2 Lattice Simulations of Vortex Formation 

We solve the system of two coupled partial differential equations ([4.2|). As in sub- 
section we find that the gradient of tachyon field becomes singular near the 



core of the defect in a finite time, forcing us to halt our lattice evolution. As in 
the case of the kink we choose as initial data T(t = 0,x,y) = Ti{x,y) = and 
T(t = 0,x,y) = Ti{x,y) given by a truncated Fourier series with random Gaussian 
coefficients with the overall amplitude small compared to b. We take H = for our 
examples, since the Hubble damping plays little qualitative role in the dynamics. In 
figures |, ||, |^ and || we plot — |T| = — ^/Tf + T2 against {x,y} with the T axis 
measured in units of b. Figure ^ shows typical initial conditions used for our numer- 
ical analysis. In figures ||, |^ and | we plot the final configurations close to t = tc 
when the gradients become infinite for various radii of compactification. (Note that 
because we are plotting — |T| rather than +|T| the vortices appear as spikes in the 
final configuration.) Again we find that vortices do form for radii of compactifica- 
tion as small as a few times M^^, even though the field is initially in causal contact 
throughout the compact space. 



4.3 Analytical Study of Vortex Formation 

As in subsection |3.3.1| we expect that in the case of the vortex there will exist fixed 
points where the field stays pinned at T = throughout the evolution. To see 
this we consider the equation of motion for Ti at t = for initial data such that 
T{t = 0,x,y) = 0. To simplify the expression we write the equation evaluated at a 
point {xq, yo) such that T{t = 0, Xq, yo) = 0: 

-Tl + din + [ 1 + (9,Tl)2 + {dy^f + {d^T^f + [dyT^f 

{d^Tif din + {dyT^f din + 2 d^n dyTi d^dyT^ 

d^n d,T2 din + dyn dyn din + {d^n dyn + dyn d^n) d.dyn = o. 
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2xR 



Figure 5: Typical initial configuration 
for numerical studies of vortex formation. 




Figure 6: The final configuration of the 
complex tachyon field for R = IM^^ . 




The equation for T2 is identical with 1^2. It is clear, then, that any point (xq, yo) 

where T/(0,xo,i/o) = 5x^/|(o,a;o,yo) = ^y^/l(o,xo,j/o)_.= ^x5s/T/|(o,xo,3^o) = (fo^ ^ = 
1,2) will be a fixed point of the evolution where Tj(t,xo,yo) = Ti(t,XQ,yQ) = for 
all t and hence the field T stays pinned at zero throughout the evolution. In the 
neighborhood of the point (xo,?/o) we thus should be able to write the field in the 
form T = u{t){x — xq) + v{t){y — yo) with u and v complex. Therefore to study 
analytically the dynamics of the field near the core of the vortex we make an ansatz 
of the type: 

Ti{t,x,y) = u{t) {x - xo), T2{t,x,y) = u{t) {y - yo). (4.3) 

We have chosen v{t) = iu{t) and u{t) real since the vortex solution is expected to 
take the form 

T (t, z,^)=u it) \[{z- Zi) J] (z - Zi) (4.4) 

i j 

where we have defined complex coordinates z = x + iy, z = x — iy. The profile 
with u{t) = const — 00 was used for the multi- vortex solutions of [2^| where 



it was shown that each holomorphic zero of (4^) corresponds to a brane and each 
anti-holomorphic zero of ([4 .41) corresponds to an antibrane. We insert now the ansatz 
( [4 .31 ) into the equations of motion ( [4 .21 ) which corresponds to studying the dynamics 
close to the core of a single vortex located at Zq = Xq + iz/Q. In this regime the 
equations of motion for the real and imaginary parts of the field, Ti and T2, give the 
same equation for the slope near z = zq: 

■■ 2 9\ 3 ?x w / , ^\ 

where we have taken H = for simplicity. We have verified numerically that this 
equation yields solutions which diverge in a finite time for generic initial data. In the 



regime where u{t) and its derivatives are large the dominant contribution to ( [4.5|) is 



2u^ X 
0^ 2 u 

which has the solution 

where the critical time when the slope diverges, tc, is fixed by the initial data. This 
singular behavior corresponds to the finite-time formation of a codimension-two brane 
in the annihilation of a brane- antibrane pair. 

4.4 Alternative Complex Tachyon Action 

The action (O) used above has the advantage of making the analysis tractable, and 



the resulting dynamics is analogous to kink formation. However, as discussed in 
subsection |4.1| this action has theoretical drawbacks. Here we consider an alternative 
description of the complex tachyon dynamics. 

The tachyon action has been calculation in boundary string field theory (BSFT) 
28| by assuming a linear tachyon profile. For a linear profile, gauge and space-time 



m 



transformations allow one to write T = UiX + i U2y, and the action one derives is 

S = -2Tp j dP+^x exp [-2W {{uixf + iu2yY)] F(4W\?) F(4W\^) (4.8) 
where the function F{z) is given by 

A'zT(z)^ 



F{z) 



2T{2z) 



To make our analysis tractable we generalize the action ([4.8|) to nonlinear profiles 
using two simplifications. The first is to consider a generalization which reduces to 
( ^4.8| ) for linear tachyon profiles only when ui = U2 = u. We feel this is a reasonable 
simplification since for the profile T = Uix + iu2y to minimize the action one requires 
both Ui and U2 to be infinite. 
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Our next simplification is to replace the complicated function F{z) by Vl + ttz. 
As justification we note that these two functions have identical behavior at large z 
since F{z) — > ^/7^z for z — >■ oo. 

We consider, then, the action 



'^Tp j y/^dP+^xexp {-\T\yb^) (l + df^Td^T) (4.9) 



^ = - 



where we have performed a rescaling of T and, for consistency with ( ^4.8| ), b = 
a/ ira'. This action has also been studied in connection with tachyon condensation 
in Writing the tachyon field in real and imaginary components as T{t,x,y) = 



Ti{t, x,y) + i T2{t, X, y) the equation of motion one derives from ( [4.9|) is 



g^'S/^VuTi + i [Ti (1 + q^'V^J^kVuTk) - 2Tk9^''V^TkV.Ti] = 0. (4.10) 
We have solved the system of two coupled partial differential equations ( [4.10| ) 



on a lattice, as in subsection |4.2| . The results are similar to those following from 
the action ( [4.1| ), though in this case the nonlinear effects are less dramatic. One 
still has defect formation, though on a longer time scale and defects begin to form 
for somewhat larger radius of compactification {R ~ 30M^^ is sufficient). The 
qualitative result that cosmic strings can form even for compactification scales well 
below the Hubble scale is unchanged. We consider now the dynamics near the core 
of the defect by plugging the ansatz 

Ti{t,x,y) = u{t){x - xo), T2{t,x,y) = u{t){y - yo) 

into ( [4.10| ) and working only to leading order in {x — xq), {y — yo). For H = the 
equation for the slope of the field at the core of the kink is 

1 

u-—u = 4.11 

leading to an exponentially increasing slope. We see that in this case the slope does 
not become singular in finite time, which can result in a different final density of 
defects. 



5. Comparison with Ordinary Cosmic String Formation 

If the potential for the tachyon has the usual runaway form, exp (— T^/6^), then once 
the field starts rolling, it continues rolling towards T = ±oo. The gradient force is 
insufficient to stop or reverse the rolling. In this section we remind the reader how 
this differs from the mechanism for formation of defects in conventional field theories, 
where their production is much more strongly suppressed. The conventional case 
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Figure 9: Comparison of the tachyon and 0^ potentials. In the case of the (p'^ theory, the 
finiteness of the slope of the kink, as well as large oscillations of the field, strongly suppress 
the formation of defects. 

corresponds to a theory with a global U (1) symmetry, the potential j [\(f)f — a^) 
and a standard kinetic term. 

In the case of 0^ theory, the evolution is well-behaved and one can follow the 
formation and the annihilation of kinks and anti-kinks indefinitely into the future. 
Comparing the two potentials exp (— and | (|0|^ — cr^) (see figure |^), we see 
that in the (f)'^ theory one expects oscillations of the field which can restore symmetry 
and wipe out the defects. 

In that case the final density of defects formed depends strongly on how fast the 
oscillations are damped, either through Hubble expansion, or through the coupling 
of the field (p to other fields. In figures ^ and |TT] we show the effect of the damping 
in the (j)'^ theory for radius of compactification R = 60M~^. A large value of H 
results not only in a higher density of defects, but also slows down the motion of the 




Figure 10: The effect of damping on Figure 11: The effect of damping on 

the formation of vortices in the (j)^ theory. the formation of vortices in the cp^ theory. 

Small damping results in a large number An unphysically large damping is used in 

of oscillations of the field, and effective order to show that the density of defects 

homogenization. which survive is larger. 
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defects that have formed, effectively reducing the rate at which they annihilate. This 
is in contrast to the tachyon field theory, in which case the Hubble damping plays 
essentially no role in determining the final density of defects. 



6. Cosmological Consequences 
6.1 String defects 

Our studies of defect formation in compact space can easily be extrapolated to im- 
ply a density of roughly one defect per string volume in the non-compact directions 
aswell. We now turn to the cosmological implications of having a very large initial 
density of strings. In normal cosmic string networks, the details of the initial con- 
ditions are not important because the network quickly reaches a scaling solution. 
This can be demonstrated quite simply, using the "one-scale" model for the energy 



density p in long strings whose characteristic length is L [p3]: 



p = -2ifp-/(P)£. (6.1) 

The terms on the r.h.s. represent respectively the effects of dilution through expansion 
and the loss due to breaking off small loops, which eventually disappear by shrinking 
and emitting gravity waves. f{P) is a function of the intercommutation probability, 
which is believed to go like / ~ a/P [^, [S^. Taking p = /i/L^ (where p is the string 



tension) and H = j3/t, one can verify that this has a stable attractor solution 

known as the scaling solution, since the string length becomes a constant fraction of 
the horizon size H~^, and the energy density in long strings also tracks that of the 
dominant component, 

/^ = :^- (6-3) 

If the initial energy density was much greater than the scaling value, we can find the 
time scale for reaching scaling by solving ( |6.1| ) in the approximation that the Hubble 
expansion term is negligible compared to the loop-emission term, giving the solution 
7(t) = 7o + 1 /(-P) ln(t/to)- Inverting, we find that the time required to reach a value 
7 = //(2 — 2f3), starting from high densities where 70 ^ 1, is 

t ^ to 6^/(1-'^) (6.4) 

which is not much greater than to- For a radiation dominated universe, with (3 = 1/2, 
the scaling solution is reached in e^/2 Hubble times in the usual case, and subsequent 
evolution is quite insensitive to the initial conditions. This conclusion is unchanged 
even for very small intercommutation probabilities. 



- 16 - 



We have investigated the approach to the scahng solution using a more detailed 
model of network evolution, which takes into account the possibility that loops may 
reconnect to long strings when the initial density is very high, and thereby possibly 
delay the onset of scaling. Suppose that the density of small loops with characteristic 
size I is defining the average distance between loops at a given time. One can 

estimate that the rate per unit volume for loops to recombine with long strings is 

^ = fx-'L-'lL^-^^. (6.5) 
at mm(x, L) 

Here / is the probability of a reconnection, and are the number densities 
of loops and long strings, respectively, frei is the relative velocity between loops and 
strings, which we take to be 0(1), and min(x,L) is the distance a loop typically 
travels before reaching a string. The probability of a collision must be proportional 
to /, not since the size of the loop in the direction parallel to the long string does 
not affect the cross section. 

In this model, long strings and small loops are treated as two separate compo- 
nents, Ps = fi/L"^ and pi = pl/x^, whose energy densities are governed by 

Ps = -2Hps - f- + pi 



p, = -3Hp, + f^-pl^- TG/z^l (6.6) 



The last term represents the power emitted by loops due to gravitational radiation, 
YGp"^ g, where T = 50 and 10"^^ < Gp < 10"^ [0-0. The loop size is taken 



to always be a fixed fraction of the long string correlation length: l/L = TGp, so 
long as this is not smaller than the fundamental string length scale Is- We have 
integrated these equations numerically, together with the Friedmann equation and 
the evolution equations for energy density in visible and gravitational radiation, 
keeping the short distance cutoff Is on the size of the loops (in fact the results do not 
change noticeably if we assume the loops remain as small as this cutoff). This more 
detailed study confirms that the scaling solution is attained in only a few Hubble 
times, as can be seen from the time evolution of the fractions of the critical density 
for each component, shown in figure p!^. We note, however, that in the above analysis 
we assumed that string velocities remain of order unity; if there is significant freezing 
out of the relative string motions in the large in the large dimensions this could have 
significant impact on the approach to scaling We note also that the effect of 
friction due to particle scattering may also significantly alter this picture [03].^ 



The initial density of strings is many orders of magnitude greater than the Kib- 
ble estimate, which gave the initial correlation length L as ~ at the moment 
of formation of the network; instead, the initial energy density of the network is 



*We thank J.J. Blanco-Pillado and Carlos Martins for bringing this to our attention. 
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(Hubble units) 



Figure 12: Fraction of critical density versus time for long strings, loops, and visible 
radiation starting from initial values ^Ig = 0.4, ^li = 0.1 and ^Irad = 0.5. Not shown is the 
fraction of energy density in gravity waves produced by decay of the loops. 



comparable to the total energy density available, ~ yU^, so that L ~ smaller 
by a factor of jj,"^/"^ /Mp than the Kibble value. Since p = p/L"^, the initial density is 
greater by a factor of M^/fi than the Kibble value. However, this huge excess is so 
shortlived that it has no observable effects. For instance, the contribution from the 
early nonscaling regime to the gravity wave background is negligible due to the small 
size of the loops which are formed and subsequently radiate during this era. Follow- 
ing ref. [^, one can estimate the amplitude of these gravity waves at frequency u 



h - M>V3 - 10 (6.7) 



as being of order 



where po is the present energy density of the universe. The estimate (|6.7|) is some 
40 orders of magnitude below the sensitivity of LIGO at the frequencies of interest, 
uj ~ 100 Hz. As for the cosmic microwave background, the wavelength of density 
perturbations created during the nonscaling regime of the network is too short to be 
relevant: initially A ~ /i"^''^, which gets stretched to the scale of the present energy 
density A ~ Po^^^ ~ 0.1 mm. This length scale is also too small to be relevant for 
the formation of primordial black holes (PBH's) since the mass contained in volume 
A^ is far below that needed for cosmologically long-lived PBH's. 

6.2 Domain Walls 

Thus the fast approach to the scaling solution for three-dimensional string networks 
erases all sensitivity to the initial conditions, even though the initial density is orders 
of magnitude greater than for conventional cosmic strings, and this is true regardless 
of the reduced intercommutation probability. However, there are situations where 
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our modified understanding of tlie network's initial conditions may make a dramatic 
difference: namely, in the case of higher dimensional defects. Let us illustrate with 
the example of D5-D5 annihilation, where two of the dimensions are wrapped on an 
internal manifold with coordinates {y 1,1/2) and the remaining three dimensions span 
the euclidean space {x^,x'^,x^). The co dimension-two defects which form from the 
annihilation are D3 branes, and these can have various orientations with respect to 
the world-volume of the parent branes. The choices are exemplified by the three 
situations: 

1. extended in X3,yi,y2 directions, localized in xi,X2: looks like Dl in 3D. 

2. extended in X2,X3,y2 directions, localized in Xi,yi: looks like D2 in 3D. 

3. extended in xi,X2,X3 directions, localized in yi,y2'- looks like D3 in 3D. 

The first case looks like ordinary cosmic strings to the 3D observer since the 
defects have only one long direction among the three large ones. Their effects have 
already been discussed. Case 3 is a network of 3-branes, all of whose dimensions 
are large. Their tension will contribute to the effective 3D cosmological constant.^ 



Case 2, illustrated in fig. |I3|, is the interesting one because these appear as domain 
walls to the 3D observer, and their energy density redshifts too slowly: p ~ in 
terms of the scale factor of the large dimensions. A single domain wall of tension 
T2 = yf within our horizon would dominate the present energy density unless r] < 1 
MeV. The effective tension of a 3-brane wrapping one compact dimension of size R 
is T2 = Rt^ ~ /i^''^; hence the string scale would have to be < 1 MeV, absurdly small. 

Our conclusion differs from that of ref. 0, which speculated that no defects 
partially localized in compact directions can form. The argument was based on the 
Kibble mechanism, and thus assumed the correlation length could not be smaller than 
H^^ ~ Mp/fi, whereas the size of the compact dimensions must be much smaller, of 
order We have seen that in fact the initial correlation length is of the same 

order as the string scale, so that this argument is invalidated. 

Still, one might be skeptical as to whether defects partially localized in the 
compact directions can survive until today, since intuitively they might be able to 



find each other and annihilate very quickly in the compact space. Ref. attempts 
to address this question, and concludes that domain wall defects like those in case 2 
cannot efficiently annihilate, since they are not completely localized in the compact 



n = -{3-pii)Hn- -^n'^ (6.8) 



(nonexpanding) space. The analysis of [^| assumes that the number density of Dp- 
branes satisfies a rate equation which is nearly the same as that governing monopoles: 

r 

where p\\ is the number of dimensions of the brane spanning the large dimensions, 

and D is the total number of spacetime dimensions. This ignores the effect of self- 

intersections for reducing the de nsity of long defects, which is known to be the 
^These can be safely assumed to annihilate quickly since they are not separated in the expanding 
space {x^,x'^,x^). 
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dominant means for string networks to reach scaling (cf. eq. ( |6.1|) ). Further, it 
unreahstically assumes that the defects are parallel, so that they will annihilate 
rather than intercommute when they meet. It is therefore not immediately clear 
how far we can trust their conclusions |37 . 

On the other hand, numerical evolution 
of domain wall networks shows that self-intersections 
are not generic, and it is suggested that the 
dominant energy-loss mechanism is direct grav- 
itational radiation rather than through colli- 
sions Furthermore it is observed that the 
approach to the scaling solution is slower for 




^ domain walls than for cosmic strings ^9 



X 



2' 



Figure 13: An array of codimension- 
two D3-branes from D5-anti-D5 anni- 
hilation, partially localized in the com- 
pact dimensions, which look like strings 
in the xi-yi plane, and domain walls in 
the large Xi dimensions. 



6.3 Full String Network Simulations 

To attempt to address the issue of whether 
domain walls disappear or not, we have con- 
sidered the dynamics of D3-branes in (5-1-1) 
dimensions in the approximation of projecting 
out the dynamics in one compact direction, 
y^, and one noncompact direction, x^, to give 
an effective 3 -|- 1-dimensional system. In this 
case the D3-branes appear as one- dimensional 
objects (see figure |T3|) and the dynamics can 
be modeled by considering string evolution in 
an anisotropic space with two large and one small dimension. In this setup those 
"strings" — string-like from the {x^,x^,y^) point of view — which span the large di- 
mensions x^ and x"^ appear as domain walls to the 3D observer while the "strings" 
which span the compact dimension appear as cosmic strings in 3D. 

Following the setup of Smith and Vilenkin |^0[ we performed numerical simu- 
lations of the defect evolution in this approximation, keeping track of the extent to 
which defects preferentially spanned the compact direction — thus appearing as 
strings in the x^, x^, x^ subspace, relative to spanning the large directions, which ap- 
pear as domain walls. Figure |1^ shows the fractional energy density in wound strings 
wrapped about the three directions function of time. This particular 

run started with equal energy in each direction so that the correlation lengths, and 
hence the interaction rates, were roughly the same. We observe that the branes wrap- 
ping the large directions lose their winding energy more quickly than those wrapping 
the compact direction: this can be attributed to several factors, including the smaller 
cross-sectional area of the long strings, and the greater energy radiated away when 
two long strings annihilate. The implication is that defects which look like strings 
to the 3D observer tend to survive preferentially over those that appear as domain 
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Figure 14: Fraction of wound string density about the three anisotropic directions (two 
large and one small) starting from equal energy density in strings wrapping each direction. 
Top curve is for strings winding in the small direction. 

walls. However since this is a toy model for the actual higher-dimensional defects, 
such a conclusion awaits validation from actual domain wall simulations. 

As further evidence that defect evolution after formation tends to favor the 
defects which are localized in the large dimensions (and hence tends to favor the 
survival of cosmic string-like defects) we consider a field theory simulation of domain 
wall evolution for a real scalar field with potential ^ (0^ "~ '^^)^ ^iiid standard kinetic 
term ^ We studied the formation and evolution of domain walls in this theory in 
an anisotropic (3 -|- l)-dimensional space with two large and one compact dimension. 
Starting from a random initial profile close to the false vacuum (as in subsection 
|3.2| ) we find that the evolution of the defects after formation tends to annihilate 
domain walls which are localized in the compact direction and to favor the survival 



of domain walls which are localized in the large directions. Figure shows a plot 
of the domain wall network late in the evolution after the domain walls localized in 
the compact direction have disappeared. 

We also explored the effect of increasing the anisotropy, varying the initial dis- 
tribution of winding modes, and varying the intercommutation probability; however, 

^Though it is of more direct interest to consider the dynamics of the tachyon field theory with 



action (3.1) the finite time slope divergence prevents us from following the dynamics aft;er the 



defects have formed. 
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Figure 15: Plot of domain wall network in an anisotropic space with two large and one 
small dimensions with L = 80Af~^. This snapshot is taken late in the evolution and shows 
that domain walls which are localized in the large dimensions are preferred by the defect 
evolution. 



the system consistently evolved to favor winding about the compact direction. These 
results seem to corroborate the claim that domain wall-like defects are suppressed. 

We note that there are several factors which we have not taken into account 
which may radically alter the picture. For example, in figure II we have unreal- 



istically neglected the dynamics in two dimensions to facilitate numerical studies. 
Furthermore, in both of the examples above we have neglected the expansion of the 
large dimensions — an effect which could play an important role in the dynamics. It is 
therefore unclear whether domain walls will pose a cosmological problem in models of 
brane-antibrane inflation where the branes driving inflation wrap the compact space. 
We feel that this is a problem which merits further investigation and it is likely that 
a complete resolution of this issue will require full simulations of 3-brane dynamics 
in an anisotropic (5 + l)-dimensional spacetime with 3 expanding dimensions. While 
such simulations would be of interest, they are beyond the scope of the present work. 

6.4 Monopoles 

Finally consider an example of how monopole-like defects may be formed through 
a cascade of annihilations in D5-D5 inflation. The initial state D5 and D5 span 
the three large dimensions and wrap two compact dimensions. These may produce 
D3 and D3 which wrap the two compact dimensions and are extended in one large 
dimension; hence these defects appear string-like to the 3-dimensional observer. A 
parellel D3-D3 pair may then annihilate to produce Dl branes and antibranes which 
can span the same large dimension as the parent D3-D3, or alternatively could wrap 
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Figure 16: Fraction of wound string density about the three anisotropic directions (one 
large and two small) starting from equal energy density in strings wrapping each direction. 



the compact dimensions. Those Dl-branes which span the large dimension appear 
as cosmic string defects to the 3-dimensional observer while those which wrap the 
compact dimensions will appear as point-like (monopole) defects. If the compact 
dimensions admit nontrivial 1-cycles (like for example) then these monopoles will 
be stable. Our results indicate that in general both string-like and monopole-like 
defects should be produced in this cascade. To understand the subsequent evolution 
of such defects we consider numerically the dynamics of Dl-branes in a 3D space 
with one large dimension and two small dimensions. Figure shows the fractional 
energy density in wound strings wrapped about the three directions x^,y^,y'^ as a 
function of time. Strings wrapping the large direction appear as genuine cosmic 
strings to the 3D observer while strings wrapping small directions y^,?/^ appear as 
monopole-like defects in 3D. As in subsection |6.2| we start with equal energy in each 
direction so that the correlation lengths are roughly equal. Again we find that strings 
wrapping the large direction lose their winding energy quicker than those wrapping 
the compact directions. Physically, this suggests that monopoles are preferentially 
produced over cosmic strings in this particular cascade of annihilations. 

The question of whether such monopole-like defects will pose a cosmological 
problem depends crucially on the long-range forces between these defects and is, 
we feel, an issue which merits further investigation (see PT| for a solution of the 
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monopole problem which is independent of inflation). 



7. Summary and Conclusions 

We have studied the formation of topological defects during the decay of a nonBPS 
brane or a coincident brane-antibrane pair. The problem was treated both analyti- 
cally, by solving the equation of motion for the tachyon field at the core of the defect, 
as well as numerically, by evolving the tachyon field on a lattice. We showed that 
defects form with a correlation length proportional to the string length rather than 
the horizon. Defects localized within compact dimensions can form even if the com- 
pactification radius is as small as 7Mg^ and the tachyon dynamics is insentitive the 
the Hubble damping. Depending on the exact form of the action {e.g., Sen's version, 
or the boundary string field theory version) the slope of the field at the defect could 
either increase exponentially in time, or else diverge within a finite time, potentially 
changing the initial density of defects. We compared the evolution of the tachyon 
field to that of a scalar field in the 0^ theory and noted that the most efficient way to 
suppress the formation of defects is through symmetry restoration, caused by large 
oscillations of the scalar field. This is not possible if the potential has a runaway 
form, as for the string tachyon, which inevitably leads to the formation of a higher 
density of defects in the string case. Once the defects form the field theory descrip- 
tion is no longer adequate, so in order to analyze the annihilation of the defects 
formed one has to use a description in terms of branes and antibranes interacting in 
a compact space. 

As a result, the initial density of string defects is much greater than previous 
estimates. For strings which are genuine ID objects, we showed that the string 
network nevertheless attains scaling behavior within just a few Hubble times, so 
that there are no observable consequences of the initial high string density. Of course 
this assumes that the network is not frustrated so that scaling can indeed be 



achieved. Whether this is the case for cosmic superstrings which are bound states of 



fundamental and D-strings is an interesting open question ^3 . 

On the other hand, we argue that in models where infiation is driven by branes 
which wrap the compact manifold (for example D5-D5), domain wall-like and monopole- 
like defects are inevitably produced. The stability and subsequent evolution of such 
defects is complicated and may depend crucially on the details of the compactifica- 
tion, for example, on whether or not the compact manifold which the parent branes 
wrap admits nontrivial 1-cycles. We leave detailed studies of whether such mod- 
els are phenomenologically viable to future investigations; however, we have shown 
that the formation of defects at the end of brane-antibrane infiation is much more 
complicated and model-dependent than one might have expected. 
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APPENDIX: Initial Conditions for Defect Formation 

Here we briefly discuss the initial conditions for defect formation at the end of brane 
inflation. During the slow roll inflationary phase the tachyon behaves as an ordinary 
massive Klein-Gordon scalar fleld (provided T <^ ^7^)- We consider here for sim- 
plicity a standard fleld theory of a scalar x in 3 + 1-dimensions whose mass squared 
parameter abruptly becomes negative. This type of theory has been considered in 
detail in [|l|, 

During the de Sitter phase (before the mass parameter becomes tachyonic) 
vacuum fluctuations yield a blackbody spectrum of produced particles {0\Nk\0) = 
^gUJk/T _ ^j^]-^ temperature T = H/{2'k) (see, for example, |^4[ for a review). 
However, in any region of de Sitter space which is small compared to the Hubble 
scale, the space is locally Minkowski, and even in the vacuum state there are quan- 
tum fluctuations quantifled by the two point functions of the flelds 

{X\k)x(k')) = ^{2iTf5\k-k') 
2\k\ 

where H = x.^ The initial stages of the string tachyon condensation are identical 
to the tachyonic preheating scenario described in [|T6], The tachyonic instability 
amplifles exponentially those modes with \k\ < m where the x fi^ld has mass squared 



parameter —m and the variance of the fluctuations grows as |T7 

{x\t)) = (x'(o)) + - dk'^;pr]^2 sinh' {tV^^F^) 

where (x^(0)) is a divergent vacuum contribution. The above result was derived in 
(3+1) dimensions, but the generalization to higher dimensions must have the form 

where rrii are the masses of the Kaluza-Klein excitations. These fluctuations grow to 
be of order the classical VEV, {x^it)) - (x^(0)) ~ m^/A before the linear treatment 
breaks down. Notice that this growth occurs on a microscopic time scale. Since in the 



Similar initial conditions are taken for defect formation at the end of inflation in 
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case of the string theory tachyon the potential has no minimum, we can conservatively 
take the initial exponential growth to be over a shorter time scale, {x^) ~ ~ M^. 
These fluctuations are much larger than the de Sitter fluctuations (x^) ~ H"^ and are 
thus the dominant seeds for defect formation. Furthermore, these fluctuations have a 
minimum wavelength comparable to the string scale which is sufficient to initiate the 
formation of defects localized in the compact dimensions. This justifies our choice 
of random initial conditions for the tachyon field in the numerical studies of defect 
formation. 
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